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Branching ratio study of ZH — > qqcc/qqbb 

Hiroaki Ono 1 * 

1 - Nippon Dental University School of Life Dentistry at Niigata 
1-8 Hamaura-cho chuo-ku Niigata, Niigata - Japan 

Precise measurement of the Higgs boson properties is an important issue of the In- 
ternational Linear Collider (ILC) experiment to verify the particles mass generation 
mechanism that the coupling strength between the Higgs boson and the fermions or 
vector bosons are proportional to the mass of each particle. Thus the measurement 
of the branching ratio of the Higgs boson an important issue to understand the mass 
of each particle. In this analysis, measurement accuracy of the Higgs boson branching 
ratio in the ZH — > qqH hadronic decay mode was studied with the cut-based analysis 
in Higgs mass of Mh ~ 120 GeV at the center-of-mass energy of ^/s — 250 GeV with 
the ILD detector model. From the analysis, we estimate the measurement accuracy of 
the relative Higgs boson branching ratio of BR(H — > cc) to BR(H — > bb) as 13.68%. 

1 Introduction 

International Linear Collider (ILC) [I] is a future 
e + e~ collider experiment for the precise measure- 
ment and the validation of the Standard Model (SM) 
physics, especially for the measurement of the Higgs 
boson property, even the discovery of the Higgs bo- 
son will be realized in Large Hadron Collider (LHC) 
experiment. In the SM, light Higgs boson mass (Mh) 
is predicted the range of Mh > 114.4 GeV from the 
study in Large Electron Positron Collider (LEP II) [2] 
and recently Tevatron experiments exclude the Higgs 
mass range around 160 < Mh < 170 GeV with the 
95% confidence level [3|. From these results, Higgs 
mass is indicated to be light (Mh < 140 GeV) and 
in this region, Higgs mainly decays to bb pair which 
forms multi-jet final state, as shown in Fig. [TJ Since 
hadron collider experiments have large QCD multi- 
jet backgrounds, the measurement of the light mass 
Higgs will not prefer in terms of the signal to noise 
ratio, ILC experiment has an advantage for precise 
measurement of light Higgs boson with large signal 

yield of multi-jet final state in lower background environment. Therefore, precise mea- 
surement of the light mass Higgs (Mh < 140 GeV) will be a primary target of the ILC 
experiment. At the lower center-of-mass energy (^/s) range around the production threshold, 
as shown in Fig [21 SM Higgs boson is mainly produced through the Higgs-strahlung (ZH) 
process, which associated with the Z boson as shown in Fig. [31 
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Figure 1: Branching fraction of the 
Higgs boson decay at the lower mass 
range (100 < M H < 200 GeV). 
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Figure 2: 
Higgs with 



Production cross section of the 




Figure 3: Higgs boson production via Higgs- 
strahlung (ZH) process. 



The largest production cross-section via 
ZH mode is obtained around the center- 
of-mass energy (y/s) at the ZH produc- 
tion threshold region as shown in Fig. [2] (a). 
Since the Z boson mainly decays to qq pair, 
the largest Higgs boson production cross- 
section via ZH process is obtained through 
the ZH —$■ qqH process. Therefore, we 
study the Higgs boson property with the 
largest production cross-section process of 
ZH — > qqH. Since Higgs boson mainly 
decays to bb pair at the Higgs mass below 
140 GeV region as shown in Fig. [2] (b), the 
final state of the ZH — > qqH process forms 
the four-jet. In ILC experiment, two de- 
tector concepts, ILD and SiD submit their 
Letter of Intent (LOI) and validated by ILC 
Detector Advisory Group (ID AG). In or- 
der to achieve the best jet energy resolu- 
tion in multi-jet environment, ILD adopt 
the Particle Flow Algorithm (PFA) suited 
detector design, which has fine-segmented 
calorimeter with strong magnetic field. In 
this analysis, we study the measurement ac- 
curacy of the branching ratio of Higgs bo- 
son with the full detector simulation for 
ZH — > qqH hadronic mode with the ILD 
detector model. 



2 Analysis tools and MC sam- 
ples 



Tabic 1: Signal and background data samples. 



Signal and Bkg 


DST data samples 


Higgs sample 


qqH, vvH, UH 


SM background 


qqqq, v£qq, vvqq, qq, ££££, gg 



For full detector simulation study, 
we use the ILD detector model 
based Monte Carlo (MC) full simula- 
tion package called Mokka, which is 
based on the MC simulation package 
Geant4 [4]. Generated MC hits are re- 
constructed and smeared in the recon- 
struction package called MarineReco 

which includes the PFA package called PandoraPFA [5]. Reconstructed his are skimmed 
and saved in the ILC common data format called LCIO. For the event analysis, we use the 
useful analysis package library called Anlib and each analysis procedure is handled through 
JSF [7] based on Root [6]. In this analysis, we assume the center-of-mass energy around the 
ZH production threshold of y/s = 250 GeV and the light Higgs mass of Mh = 120 GeV. 
Each data sample is scaled to the integrated luminosity of C = 250 fb _1 and the beam 
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polarization to P(e + ,e ) = (+30%,— 80%). The main backgrounds for ZH hadronic de- 
cay mode are considered as following processes: ZH — > Z* /^ — > qq, e + e~ — > WW/ZZ — > 



qq'q"q'" or qqq'q' 



iiqq, vvqq, 



WW -> vlqq' and e+e~ -> ZZ 



Visible energy 



■ H->cc 

— ZH^qqh 

— ZH-»llh 

— ZH-»nnh 




Figure 4: Visible energy distribution to select 
the Z — > qq hadronic mode from pre-mixcd 
sample. 



3 Event reconstruction 

Since the final state of the ZH —$■ qqH 
mode forms four-jet, after the PandoraPFA 
clustering, forced four-jet clustering based 
on Durham jet-clustering algorithm has ap- 
plied. In order to select the best jet 
pair combination from the four-jet, d — 
(M 12 - M z f + (Af 34 - M H f value is eval- 
uated, where M12 is a Z candidate di-jet 
mass, M34 is a Higgs candidate di-jet mass, 
AIz/h ar e Z ancl Higgs boson masses. From 
the four-jet, minimum d value jets combi- 
nation is selected as best Z and H candi- 
dates. In order to improve the background 
reduction from WW/ZZ with mass distri- 
bution, kinematic fitting is applied with 
following constraints after the jet pairing; 
£ Ei - Ecu = 0, £ Pi = 0, |Mia - M 34 | = 
\Mz — Mh\ (di-jet mass difference consis- 
tent with Z and H mass difference), where 
Ei and Pi are i-th jet energy and momen- 
tum which is sorted by energy, respectively. 

After the jet pair combination and kinematic fitting, event selections are applied. For 
the ZH hadronic decay mode sample (ZH — > qqH) selection from the pre-mixed sample 
of ZH — > qq/i/p/liH, following signal classification is applied: (0). Visible energy cut : 
E V i S > 170 GeV and no high momentum tracks (Pi ep ton > 15 GeV), After the classifica- 
tion, we apply the following selection criteria. In order to select the four jet reconstructed 
events, we require the number of charged tracks in each jet above four (N c h a rged > 4) and 
logarithm of the jet reconstruction Y-value threshold from three to four jet (Y34) should be 
— logY34 > 2.7. After that, thrust, thrust angle [cosOthrust) and jets angle (0h,z) cuts are 
required to suppress the ZZ background from the difference of sphericity of final state jet 
shape. Finally we require the consistency of mass of the reconstructed jets pair after the 
kinematic fit that Z candidate jets pair should be consistent with Mz and the other pair 
should be consistent with Mh ■ For the reduction of hard photons and ISR photons, finally 
highest photon energy cut is applied. Table. [5] shows the background reduction summary 
after the selections. 

4 Branching ratio measurement 

Branching ratio (BR) of Higgs boson is related to the mass of the fermions and gauge bosons. 
After the background reduction, we evaluate the measurement accuracy of the branching 
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Table 2: Background reduction summary and its efficiency. 




Selection criteria 


H — »• cc 


H -^bb 


Higgs Bkg 


SM Bkg 


Before the classification 


2914 


53480 


23447 


53333000 


After the classification 


1693 


29075 


9198 


20528900 


^ charged -^ ** 


f238 


22204 


5721 


3323060 


- log(V34) > 2.7 


1218 


21869 


5694 


2635920 


thrust < 0.95 


1217 


21858 


5693 


2584510 


\cos6 t hrust\ < 0.96 


1157 


20831 


5427 


2295690 


f05° < 6 H < 165° 


1080 


19393 


4941 


1908300 


70° <6 Z < 160° 


1028 


18490 


4705 


1776150 


f fO < M HfH < 140 GcV 


982 


17666 


4411 


1209100 


80 < M Zfu < HO GcV 


982 


17665 


4409 


1206570 


E 1 < 20 GcV 


895 


16288 


4063 


1036990 


Efficiency after classification 


52.9% (e cc ) 


56.0% (£ 66 ) 


44.2% 


5.0% 



ratios (BR). In this analysis, we evaluate the relative branching ratios of H — > cc to H — > bb: 

-=-. In order to evaluate the measurement accuracy of BR, we apply the flavor- 

BR[H — > bb) 

likeness template fitting [9] after the all event selections. For each flavor, neural-net training 

is performed using Z — »• qq Z-pole (y^s = 91.2 GeV) samples in LCFIVcrtex package. Here, 

6c-likeness is a c-likeness trained only with Z —> bb background. Each flavor-likeness for 

di-jct is defined as: x- likeness = ; — -, where x± 2 arc the neural- net 

X\X2 + (1 — #l)(l — X2) 

trained output for 6, c and be flavor in each jet from the vertcxing package of LCFIVcrtex [8] 
in ilesoft. The analysis procedure of the template fitting is as following: (1). Prepare the 
6, c and be three dimensional flavor-likeness template samples for H — > bb, cc, others and 
SM Bkg decay modes; (2). Apply the Toy-MC template fitting test to evaluate the r^, 
r C c, fathers and r^kg, where r^, r cc are the number of entry ratios of H — > bb, cc after the 
selection cuts to the entry predicted from SM Higgs branching ratios, and r ot hers is a ratio 
of another Higgs decay modes which is described as r ot hers = 1 — n* — fee, and r^ g is a 
normalized factor for the entries of SM background. (3). Estimate the relative branching 
fraction of BR(H — > bb)/BR(H — > cc). In order to evaluate the fractions of Higgs decay 
mode in each sample, template fitting is applied with minimizing following \ 2 value: 



N data 



EEE 

i=\ j=l fe=l 



E 

-bb.ee. others 



r s 



N 



ZH 



N> 



■N 



1 Jk - r bkg ■ n:, 



bkg 
ijk 



Ajall 
ijk 



(1) 



N$ 



arc 



where r s represents the fitted parameters of r^, r cc , r others and r^g. iV?- fe 

the number of expected entries in each 3D sample bin (i, j, k), since each histogram is 

separated into rn,, n c and rif, c bins. Nf°£ a are the number of simulated entries in each bin 

(i, j, k) by Toy-MC, which fluctuated with the Poisson distribution in each samples. 

In order to estimate the rbb and r cc , template fitting is applied with fluctuating the data 

by Poisson distribution and perform the 1,000 times Toy-MC analysis. Fitted rbb and r cc 
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Figure 5: Image of the two-dimensional template sample histograms in b- likeness and c- 
likeness. 



fractions are obtained from the distribution of the Toy-MC template fitted with Gaussian, 
as shown in Fig [6] 



From the 1,000 times Toy-MC 
template fitting, we obtain the Tbb 
and r cc as 0.767 ± 0.002, 0.422 ± 
0.006, which reproduce the true 
signal fractions of r^b (0.765) and 
r cc (0.0422), respectively. 
From the limitation of the com- 
monly reconstructed full simula- 
tion samples, especially for SM 
background, number of entries 
in each template sample 3D his- 
togram, which depends on the his- r . n j ,. , ., c_ ,, , , 

' . . , r igurc o: Tbb and r cc distribution from tnc template 

togram binning, become a system- „,,. „ xrt ~. 
? . . , n i httmg loy-MC. 

atic uncertainty of the fitted pa- 
rameters. In order to reduce the 
histogram binning dependence, we 

apply the histogram smoothing for template samples. Fig. [7] shows the binning dependence 
after the smoothing. 
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5 Results 



The relative branching ratio in H — > cc to H 
evaluated with following equation: 

BR(H -> cc) r cc /e cc 



BR(H -> 66) r 6b /e: b6 ' 



66 is 



(2) 



where the £bb/cc arc the efficiency of H — > bb/cc events 
after the selections as shown in Tabled] From the Eq. @, 
relative branching fraction is obtained as: 

— ZtS^L = 0.058 ± 0.008, which corresponds to the 

BR{H -> 66) ' P 

measurement accuracy of 13.68%. 

6 Conclusion 

Measurement accuracy of Higgs branching ratio in Higgs 
hadronic decay mode ZH — > qqH in ILC experiment is 
evaluated at the y/s = 250 GeV with assuming the Higgs 
mass of Mjj = 120 GeV and the integrated luminosity 
of C = 250 fb _1 . From the template fitting analysis, 
measurement accuracy of the relative branching ratio of 
H — > cc to H — >• 66 is evaluated as 13.68%. 
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Figure 7: Binning dependence of 
the uncertainty of r&(,/ cc and r\,kg 
after the smoothing. 
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